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Abstract 

A next-to- leading order QCD analysis of spin asymmetries and structure functions in polarized 
deep inelastic lepton nucleon scattering is presented within the framework of the radiative 
parton model. The Q 2 -dependence of the spin asymmetry A\{x, Q 2 ) is shown to be non- 
negligible for x-values relevant for the analysis of the present data. In the second part of the 
paper we present LO results for radiatively generated spin-dependent parton distributions of 
the photon and for its structure function gj. 
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1 Introduction 



The past few years have seen much progress in our knowledge about the nucleons' spin struc- 
ture due to the experimental study of the spin asymmetries Af(x, Q 2 ) ~ gi(x, Q 2 ) / 'F± (x, Q 2 ) 
(N = p, n, d) in deep-inelastic scattering (DIS) with longitudinally polarized lepton beams 
and nucleon targets. Previous data on A\ by the SLAC-Yale collaboration |1| have been 
succeeded by more accurate data from [2-4], which also cover a wider range in (x, Q 2 ), and 
results on A™ and Af have been published in || and || |7|, respectively. 

On the theoretical side, it has become possible to perform a complete and consistent study 
of polarized DIS in next-to-leading order (NLO) of QCD, since the calculation of the spin- 
dependent two-loop anomalous dimensions, needed for the NLO evolution of polarized parton 
distributions, has been completed recently ||. A first such study has been presented in M, 
where the underlying concept has been the radiative generation of parton distributions from a 
low resolution scale /i, which in the unpolarized case had previously led [|H|] to the remarkably 
successful prediction of the small- x rise of the proton structure function F$ as observed at 
HERA |fTl]| . The main findings of this NLO analysis |§, which followed the lines of an earlier 
leading order (LO) study |I2]| , will be collected in section 2. 

New precise data on polarized DIS will be added in the near future from the HERMES 



experiment |13| at HERA. Moreover, it is no longer inconceivable to longitudinally polarize 
HERA's high-energy proton beam. The corresponding situation with unpolarized beams has 



already demonstrated [|11]] that at such high energies the ep interactions, since dominated by 
the exchange of almost real (Weizsacker- Williams) photons, can reveal also information on 
the parton content of the photon in addition to that of the proton. Therefore, a polarized 
ep-collider mode of HERA could in principle serve to explore the spin-dependent parton 
distributions of circularly polarized photons. These are completely unmeasured and thus 
unknown up to now, and one has to invoke models 0] in order to study their expected 
size and to estimate the theoretical uncertainty in predictions for them. For this purpose, 
it seems worthwhile to also resort to a radiative generation of the photon's polarized parton 
distributions |TJ|], since the corresponding predictions for the unpolarized photon El have 
again been phenomenologically successful [jnj]. This topic will be covered in section 3. 
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2 NLO Radiative Parton Model Analysis of Polarized 
DIS 



Measurements of polarized deep inelastic lepton nucleon scattering yield direct information 
[1-7] on the spin-asymmetry 

aN( q2\ ^ 9i( x , Q 2 ) = 9i{x,Q 2 ) /.x 

1 ' V ' ~ F?(x, Q 2 ) F 2 N (x, Q 2 )/ [2x(l + R N (x, Q 2 ))] ' 1 } 

where N = p,n,d and R = F L /2xF 1 = (F 2 - 2xF 1 )/2xF 1 . In NLO, gf {x,Q 2 ) is related to 
the polarized (Af N ) quark and gluon distributions in the following way: 

g?(x,Q 2 ) = lT,<\^ N (x,Q 2 ) + Aq N (x,Q 2 ) + 



+ 



2 , 

a s (Q' 
2tt 



AC q *(Aq N + Aq N )+AC g *Ag \ (2) 



with the convolutions (*) being defined as usual, and where the appropriate spin-dependent 
Wilson coefficients AC; in the MS scheme are given, e.g., in ||. The NLO form of the 
unpolarized (spin-averaged) structure function F^(x,Q 2 ) is similar to the one in (2) with 
Af N (x,Q 2 ) — > f N (x,Q 2 ) and the unpolarized Wilson coefficients given, for example, in flT6f . 

The NLO (Revolution of the spin- dependent parton distributions Af(x, Q 2 ) = Af p (x, Q 2 ) 
is performed most conveniently in Mellin n-moment space where the solutions of the evolution 
equations (see, e.g., refs. [16-18]) can be obtained analytically, once the boundary conditions 
at some Q 2 = fi 2 , i.e. the input densities Af(x, /i 2 ) to be discussed below, are specified. These 
(Revolutions are governed by the spin-dependent LO JHJ and NLO (MS) anomalous di- 
mensions. Having obtained the analytic NLO solutions for the moments of parton densities 
it is simple to (numerically) Mellin-invert them to Bj0rken-x space as described, for example, 
in |17|, As seen in (2), the so obtained Af(x,Q 2 ) are then convoluted with the Wilson 
coefficients ACi to yield the desired g±(x, Q 2 ). 

To fix the polarized NLO input parton distributions Af(x,Q 2 = fi 2 ) we perform fits to 
the directly measured asymmetry A±(x, Q 2 ) in (1), rather than to the derived gi(x,Q 2 ), 
mainly because for the experimental extraction of the latter often the assumption of the 
Q 2 -independence of A±(x, Q 2 ) is made, which is theoretically not justified as we will see 
below. As mentioned in the introduction, the other main ingredient of our NLO analysis 
@ is that we follow the radiative (dynamical) concept |I| by choosing the low input scale 
Q 2 = fx 2 = 0.34 GeV 2 and implementing the fundamental positivity constraints 

\Af( Xl Q 2 )\<f(x,Q 2 ) (3) 

down to Q 2 = fi 2 . Therefore we shall use all presently available data [2-7] in the small-x 
region where Q 2 ^ 1 GeV 2 , without introducing lower cuts in Q 2 as was usually necessary in 
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previous analyses PC| ]. A further advantage of this approach is the possibility to study the 
Q 2 -dependence of A±(x, Q 2 ) over a wide range of Q 2 which might be also relevant for possibly 
forthcoming polarization experiments at HERA. The analysis affords some well established 
set of unpolarized NLO parton distributions f(x,Q 2 ) for calculating Ff(x,Q 2 ) in (1) which 
will be adopted from ref . HlO| . 

In addition to (3), the polarized NLO parton distributions Af(x,Q 2 ) are, for the SU(3)f 
symmetric 'standard' scenario, constrained by the sum rules 

1 dx (Au + Aw - Ad - Ad) (x, f i 2 )=g A = F + D = 1.2573 ± 0.0028 (4) 
J dx (Au + Au + Ad + A d - 2(As + As)) (x, /i 2 ) = 3F - D = 0.579 ± 0.025 (5) 



with the values of and 3F — D taken from ||21| . It should be noted that the first moments 
of the flavor non-singlet combinations which appear in (4) and (5) are Q 2 -independent also 
in NLO H, |,[23fl. 

As a plausible alternative to the full SU(3) f symmetry between charged weak and neutral 
axial currents required for deriving the 'standard' constraints (4) and (5), we consider a 
'valence' scenario [12, ^] where this flavor symmetry is broken and which is based on the 
assumption |24| that the flavor changing hyperon /3-decay data fix only the total helicity of 
valence quarks: 

[ dx(Au v - Ad v )(x,fi 2 ) = g A = F + D = 1.2573 ± 0.0028 (4') 
Jo 

I dx(Au v + Ad v )(x,fi 2 ) = 3F-D = 0.579 ±0.025 . (5') 
Jo 

We note that in both above scenarios the Bj0rken sum rule manifestly holds due to the 
constraints (4), (4'). Our optimal NLO input distributions at Q 2 = jj 2 = 0.34 GeV 2 subject 
to these constraints can be found in ||. 

A comparison of our results with the data on Af(x, Q 2 ) is presented in Fig.l. Obviously, 
very similar results are obtained for the two scenarios considered above. The Q 2 -dependence 
of Ai(x, Q 2 ) is presented in Fig.2 for some typical fixed x values for 1 < Q 2 < 20GeV 2 relevant 
for present experiments. In the (x,Q 2 ) region of present data [2-7], Ai(x,Q 2 ) increases with 
Q 2 for x > 0.01. Therefore, since most present data in the small- a: region correspond to small 
values of Q 2 ^ 1 GeV 2 , the determination of g\ (x, Q 2 ) at a larger fixed Q 2 (5 or 10 GeV 2 , 
say) by assuming A±(x, Q 2 ) to be independent of Q 2 , as is commonly done [2-7], is misleading 
and might lead to an underestimate of g\ by as much as about 20%. Results for the (also 
non- negligible) Q 2 -dependence of A\ can be found in || [12| . The assumption of approximate 
scaling for Ai(x, Q 2 ) is therefore unwarranted and, in any case, theoretically not justified 
liii- 
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Figure 1: Comparison of our NLO results for Af(x, Q 2 ) as obtained from the fitted inputs at 
Q 2 = /i 2 for the 'standard' and 'valence' scenarios with present data [2-7]. The Q 2 values adopted 
here correspond to the different values quoted in [2-7] for each data point. 
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Figure 2: The Q 2 -dependence of A\(x, Q 2 ) as predicted by the NLO QCD evolution at various 
fixed values of x. The LO results are from [12]. 
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A further result of our analysis is that the polarized gluon density Ag(x, Q 2 ) is hardly 
constrained by present experiments. Similarly agreeable fits as those shown in Fig.l to all 
present asymmetry data can also be obtained for a fully saturated gluon input Ag(x,[i 2 ) = 
g(x,fi 2 ) as well as for the less saturated Ag(x,fi 2 ) = xg(x,[i 2 ) or even a purely dynamical 
input Ag(x,fi 2 ) = 0. We compare such gluons at Q 2 = 4 GeV 2 in Fig. 3. The variation of 
Ag(x, Q 2 ) allowed by present experiments is indeed sizeable. This implies, in particular, that 
the Q 2 -evolution of gi(x,Q 2 ) below the experimentally accessible x-range is not predictable 
for the time being. 



To conclude this section, our results H [12|] demonstrate the compatibility of the restric- 
tive radiative model [n| with present measurements of deep inelastic spin asymmetries and 
structure functions. 
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Figure 3: The experimentally allowed range of polarized gluon densities at Q 2 = 4 GeV 2 for the 
'valence' scenario with differently chosen Ag(x,fj, 2 ) inputs. 



3 Spin-dependent parton distributions of the photon 

As mentioned in the introduction, the photon's polarized parton distributions, A/ 7 = Ag 7 , 
Ag 7 , are presently unmeasured and thus completely unknown. Theoretical expectations can, 
however, be derived |L4j| by assuming a radiative generation of the A/ 7 along the lines followed 



[p~5f successfully for the unpolarized photonic parton distributions, / 7 = q' y ,g' y - In [|TjJ a VMD 
valence-like structure at the low resolution scale Q 2 = [i 2 = 0.25 GeV 2 was imposed as the 
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input boundary condition, assuming that at this resolution scale the photon behaves like 
a vector meson, i.e., that its parton content is proportional to that of the p-meson. Since 
nothing is known experimentally about the latter, the parton densities of the neutral pion 
as determined in a previous study |27j were used instead which are expected not to be too 
dissimilar from those of the p. Unfortunately, this procedure is obviously impossible for 
determining the VMD input distributions A/ 7 (x, /i 2 ) for the polarized photon. Some help is 
again provided by the positivity constraints 

\Ar(x,p 2 )\<p(x,t, 2 ) , (6) 

and it is interesting to see how restrictive these general conditions already are. For this 
purpose we consider |14[ two very different scenarios with 'maximal', A/ 7 (a;, p 2 ) = / 7 (a:,p 2 ), 
and 'minimal', A/ 7 (a;, p 2 ) = 0, saturation of (6). We mention that a sum rule expressing the 
vanishing of the first moment of the polarized photon structure function gj was derived from 
current conservation in [28] which, in the LO considered here, is equivalent to 

/ 1 Ag 7 (x,p 2 )=0 . (7) 
Jo 

This sum rule can in principle serve to further restrict the range of allowed VMD inputs. 
On the other hand, we are interested only in the region of, say, x > 0.01 here, such that 
the current conservation constraints (7) could be implemented by contributions from smaller 
x. To estimate the uncertainties in the predictions for the A/ 7 (x, Q 2 ) stemming from the 
insufficiently known VMD input we therefore stick to the two extreme scenarios discussed 
above, even though strictly speaking the maximally saturated input violates the sum rule (7). 

Starting from the two different boundary conditions for A/ 7 (x, /i 2 ) we generate A/ 7 (x, Q 2 ) 
at Q 2 > fi 2 as in |15| replacing the unpolarized splitting functions in the (Revolution equa- 
tions by their polarized counterparts ||19|| . In view of the uncertainties in the input, we restrict 
our calculations to the leading order, although in principle a NLO analysis has become possi- 
ble now pn] by exploiting the results for the polarized two-loop splitting functions of ref. [BJ. 



The resulting parton asymmetries \T% 

are shown in Fig. 4 for Q 2 = 30 GeV 2 . As can be seen, there are quite substantial differences 
between the results for the two scenarios. Fig.5 shows the corresponding LO results for the 
polarized photon structure function 

S 7 (*,Q 2 )=$> 2 Ag 7 (*,Q 2 ) , (9) 
q 

where charm contributions from the subprocesses 7*7 — > cc and 7*^ — > cc have been included 
(see |3(| for further details). For comparison we also include in Fig.5 the 'asymptotic' result 
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Figure 4: The predicted photonic parton asymmetries Aj(x, Q 2 ) at Q 2 = 30 GeV 2 as defined in 
(8) for 'maximal' (solid lines) and 'minimal' (dashed lines) saturation of the VMD input. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

X 

Figure 5: LO predictions for the photon's spin-dependent structure function xgi(x,Q 2 )/a at 
Q 2 = 10 GeV 2 for 'maximal' (solid) and 'minimal' (dashed) saturation of the VMD input. The 
dotted line shows the result for the LO 'asymptotic' solution for 3 active flavors also considered in 
[31]. 
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for gl which was also considered in refs. [31]. g\ would in principle be accessible in polarized 



e + e collisions |Kj or in erf processes. In a polarized collider mode of HERA, the spin- 
dependent parton distributions of the photon will show up in processes like jet [[32] or 



heavy flavor photoproduction when the photon is resolved into its hadronic structure [^9| . 
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